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Abstract 

1^ I Distribution functions for random variables that depend on a parameter are com- 

,J^ ' puted asymptotically for ensembles of positive Hermitian matrices. The inverse Fourier 

C^ . transform of the distribution is shown to be a Predholm determinant of a certain op- 

erator that is an analogue of a Wiener-Hopf operator. The asymptotic formula shows 
that up to the terms of order o(l), the distributions are Gaussian. 

>. 

^ 1 Introduction 

t^^ I In the theory of random matrices one is led naturally to consider the probability distribution 

^-^ ' on the set of eigenvalues of the matrices. For NxN random Hermitian matrices one can show 

O ' that under reasonable assumptions, the probability density that the eigenvalues Ai, . . . , Aat 

r-| ■ lie in the intervals 

c3 ! {xi,xi + dxi),. . . ,{xn,xn + dxN) 



is given by the formula 



> 

X' Pn{xi,...,xn) = ^7:detK{xi,Xj)\^--^ (1) 

where 

N-l 

KN{x,y) = J2Mx)My) (2) 

1=0 

and (f)i is obtained by orthonormalizing the sequence lx^e~^ /^ > over R. 

For NxN positive Hermitian matrices the probability density has the same form except 
that (pi is replaced by the functions obtained by orthonormalizing the sequence <x'^/^e~^/^x' > 
over R"*". We will not describe here exactly how these particular densities arise but instead 
refer the reader to ISll. 



* ebasor@calpoly.edu. Supported in part by NSF Grant DMS-9623278. 



We can define a random variable on the space of eigenvalues by considering f(xi, . . . ,xn) 
where / is any symmetric function of the Xj's. A particular case of interest is a random 
variable of the form YliLi fi^i)-, where / is a function of a real variable. Such a random 
variable is generally called a linear statistic. 

In previous work |]^, ^, |I|, ^ , the variance of the random variable was computed in the 
large N limit. More precisely, the function / and the kernel Ki\i{x,y) were suitably rescaled 
so that the limit as A^ — * oo of the variance could be computed. The precise details of this 
are in the next section. 

Our goal in this paper is to compute the distribution function for a class of the linear 
statistics that depend on a parameter a. We now describe the sections of the paper and 
main results. In the next section we outline the random matrix theory and show how the 
distribution functions can be computed using Fredholm determinants. In section 3 we replace 
the function f{x) in the linear statistic by /a(x) = f{x/a). For random variables of this 
type we show that the inverse Fourier transform of the distribution function 0(A;) has an 
asymptotic expansion of the form 

0(A;) ~ 6"'='+''^ (3) 

as a — i> oo. This of course implies that the actual distribution is asymptotically Gaussian. 
Here a and b depend on / and a. This is proved for both the Hermitian matrices and positive 
Hermitian matrices. In the latter case with u = —1/2, a very simple proof is given in Section 
3. For z/ > —1/2, a completely different proof is obtained in Section 4. 

Most of the results are obtained by using simple operator theory identities in the the- 
ory of Wiener-Hopf operators. The central idea is that the various quantities which yield 
information about random variables can all be computed in terms of traces or determinants 
of integral operators. Some of the computations lead directly to a familiar problem in the 
theory of Wiener-Hopf operators, while others require modifications and generalizations of 
these results. 

2 Preliminaries 

In this section we show how to compute the mean, variance, and inverse Fourier transform 
of the distribution of the random variable. Computations for the mean and variance have 
been given before in many places. However, we reproduce all of these here for completeness 
sake and also to highlight the use of operator theory ideas. 

We begin by considering P/v ioi N x N random Hermitian matrices. We want to consider 
large matrices and thus we let A^ -^ oo, but this leads to a trivial result unless we rescale 
Kn in a particular way. We replace Kn{x, y) with 




-.,■,, (4) 

'2N \V2N V2N 



Rescaling Kn is equivalent to rescaling the mean spacing of the eigenvalues. (See [0 for 
details.) 

From the theory of Hermite polynomials it is easy to see that as A^ ^ oo, 

T / X y \ sm{x-y) 



This last function is known as the sine kernel. Now consider a random variable of the form 

N 

J2fix,V2N) 

i=l 

where in all that follows / is a continuous real- valued function belonging to Li (R) and which 
vanishes at ±00. The appearance of the y/2N should not be surprising here since the above 
rescaling spreads out the eigenvalues and hence should be reflected in the random variable. 
The mean fi]y is 

N 



I '" J Yl fixiV2N)PN{xi, ..., XN)dxi ■ ■ ■ dxN- (6) 

Now the function P^ has the important property [§] 
N\ f f 

.j^ _ X, / ■ ■ ■ / ^NiXi, . . . , X„, Xn+1, ..., XN)dXn+l " " " dXN = dct K{Xi, Xj) |^^-^ . (7) 

Thus, d^) is easily seen to be 

/•oo 

f{xV2N)KN{x,x)dx (8) 



which, after changing x to x/y/2N, becomes 

1 IX V \ 

fix) , — -.Kn , , , dx. 

.^ ' V2N \V2N V2NJ 

Thus, as A^ — i> oo, 

/■oo 

/iAT — > /i = / f{x)K{x,x)dx (9) 



where K{x, y) is the sine kernel. 

A very similar computation for the variance vaiNf, again using (0), yields 

var/ := hm var„/ = - / //(x)/(y)A-(.,,)dxdy + I fix)Kix,x)dx. (10) 

Both the mean and the variance can be interpreted as traces of certain Wiener-Hopf opera- 
tors. To see this, consider the operator A{f) on L2(— 1, 1) with kernel 

1 r°° 
- f{t)e-<^-^Ut. (11) 

ZTT J-oo 

This operator can easily be seen to be the product J-'AifJ-'^^P where Pg = X(-i,i)fi') -^/fl' = 
fg and JF is the Fourier transform. A moment's thought shows that n = tr{y4(/)} and 

var/ = tr{A(r)-(A(/))n. 

A more difficult, yet also straightforward problem, is to find an expression for the distri- 
bution function of a random variable of this type. A fundamental formula from probability 
theory shows that if we call the probability distribution function 07v, then 



')N 



{k)= r ■■■ r e''^^r=i^(":'^^)p^(xi, . . .,XN)dxi ■ --dxN. (12) 

J —oo J —oo 



Thus, 

/OO fOO -" 

■ ■ ■ / n e'''^^''^'^^ Pn{x^, ...,XN)dxi--- dxN 
-OO J ~oo •_-, 

/OO /"OO " , 

■ ■ ■ / n((e*'^^"^^ ^ - 1) + l)^Jv(a:i, . . . ,a;;v) c^^i ■ ■ ■da:jv 
-OO -/— OO ■_-! 

/OO TOO ^ , -^^ , , 

. . . / {1 + y^(gifc/(:r,V2JV) _ 1) + y^(^^i.kf(x,V2N) _ ^-^^^ckfixiVm) _ ;L) _^ ^ ^ ^| 
-OO J— OO •_i • ^7 



^-oo ^.^^ ^.<; 

xP/v(a^i5 • • • , xn) dxi ■ ■ ■ dxjy 

1 + - / (e^'=/(-^^) - l)i^^(x, x) dx 

1! J — CO 



T /"OO /"OO ^ ^ 

Zl J —OO J —OO — ~ 

1 fOO noo " ^ 



In each integral we rescale to obtain 

][ ^oo ]^ /-OO roo 



I /"OO I /"OO roo 

(pNik) = 1 + TT / K'{xi,xi)dxi + — / K'{xi,X2)dxidx2 

1! J — CO 2,\ J —CO J —CO 

T /"OO /"OO 

H ^ l^i / ■ ■ ■ / ^'(^1' ■■■,XN)dxi--- dxN (13) 



OO J —OO 



where 

Letting A^ — > cxd we see this is the formula for the Fredholm determinant det(/ + K) where 
K has kernel 

i^(a;,y) = (e-^/W-l) ^^^/^~f . (15) 

7v{x - y) 

As before we can express this last quantity in terms of the operator A{a) 

(j){k) = lim (j)N{k) = det(/ + A{a)) (16) 

Af^oo 

where a{x) = e^^^^''^ — 1. 

The preceding computations can all be carried out in the case of positive Hermitian 
matrices. In this case we replace K]\j{x,y) with 

AN ^AN' an' 
and from the theory of Laguerre polynomials we see that as A^ -^ oo 

4iV^^^4iV' 4iV^ ^ 2{^^) ^^^^ 



where Jy is the Bessel function of order v. The details of this are found in |T3[. The rescahng 
here forces the eigenvalue density to be bounded near zero and is called "scaling at the hard 
edge." The kernel ( p!7D is known as the Bessel kernel. 

We can again write the mean, the variance, and the Fourier transform of the distribution 
in terms of operators. This time the relevant operator B{f) is defined on 1/2(0, 1) with kernel 
given by 

K{x,y)= t^f{t).Utx).h(ty)dt. (18) 



If we begin with the linear statistic (the y/x is merely for convenience, and we again assume 
that / is continuous, in Li(R+) and vanishes at +0x3) 



N 



E/(V^^4iV), 



(19) 



then nearly identical computations show that 

/i = iiB{f) 
var/ = tT{B{p)-{B{f)f} 
4){k) = det(/ + 5(a)) 

where a = e**'-^'-^-' — 1. We summarize these results in the following: 

Theorem 1 (a) Given a random variable of the form YliLi fixi\/2N) defined on the space 
of eigenvalues of N x N Hermitian matrices with probability distribution given in (0j, we 
have 



/i 
var/ 

m 



limjv^oo fJ'N 
limAf_oo var at/ 

limAr^oo07v(fc) 



tr(A(/)) 

tr {A{f') - {A{fm 

det(J + A{a)) 



where <j{x) = e'^^^^^'' — 1. 

(b) Given a random variable of the form Y^f=i f{^/xiAN) defined on the space of eigen- 
values of positive N x N Hermitian matrices, we have 



var/ 

m 

where aix) = e^^^^^"* — 1. 



limiV-+oo ^^N 

limAT^ooVarjv/ 
YinvN ^00 4>N{k) 



tr{B{f)) 

tr {B{f') - {Bifm 

det{I + B{a)) 



When linear statistics are considered 0, |T^, one is often concerned with a statistic of 
the form J2^i fi^i/oi) where a is a real parameter approaching infinity. This is the case, for 
example, in the study of disordered conductors where large a corresponds to a high density 
metallic regime. The above formulas still hold, of course, but now they depend on the 
parameter. We will call the operators that depend on the parameter a by Aa{f) and Ba{f), 
respectively. In the next sections we will compute the mean, variance, and distribution 
function asymptotically as a 



00. 



3 The Mean, Variance, and Distribution Function as 

a ^ oo 

For random Hermitian matrices, computing the various limits are applications of the con- 
tinuous analogues of the Strong Szego Limit Theorem. For then, Aa{f) is just the classical 
Wiener-Hopf operator defined on the interval {—a, a), and all of the quantities are known 
asymptotically as a ^ oo. We provide the answers here for completeness. 

Theorem 2 Assume that f G Li{lV) is continuous, and vanishes at ±00 and that in addition 
its Fourier transform f satisfies 



|x||/(x)|^ t/x < 00. 

-00 

Then 



zvr J-00 



and 



var/ = 2 / xf{x)f{—x)dx + o(l) 
Jo 



a rao ^ rao 



(j){k) ~ exp < — / ikf{x) dx — k / xf{x)f{—x) dx \ . 
.zvr J-00 JQ 



The Bessel case is significantly more complicated. There is no corresponding Szego type 
theorem. We begin by computing the mean. The operator Baic) has kernel 

'xytf{t/a)Ji,{tx)Jy{ty) dt. 



/o 
Thus the mean /x is given by 

nl /"OO 



H = xtf{t/a)J^{tx) dtdx 

Jo Jo 

/•OO /•! 

= a^ / xtf{t)Jl{atx)dxdt 
Jo Jo 

a" 



Now 



/ f{t) xtJl{atx)dxdt. (20) 

JO Jo 

1 1 , , 

xJl{atx) dx = -\ J^{at) — Jy^i{at)Jy^i{at) \ 



and 



2i/ 

Jy-i[at) = -J„+i{at) H -J^{at). 

at 

Therefore the integral (^0|) becomes 

2u 



a 



/•OO fyf f '2,1/ ^ 

Jo ^^^^T |'^'("^) + ^'+iM - —Ju+Mt)J,{at)^ dt 



6 



or 

at 



a / fit)— {J^iat) + J,'+i(at)} dt - av f{t)J,+i{at)J,{at) dt. (21) 



/o 2 ■- ' ' ■ J jQ 

The first integral equals 

- r fit) dt + 0(1) (22) 

n Jo 

which can be easily seen by using the asymptotic properties of Bessel functions. The second 
integral is asymptotically 

^/(0) + o(l). 

This uses the identity /q°° Ji,+i{x)Ji^{x) dx = \. 
Thus we have 

/.= - r/(t)rft-^/(0) + o(l). (23) 

TT Jo 2 

For the variance we refer to [Q] where the calculation was already done. There it was 
shown that 

var/ ~ — / \M{f){2iy)\^yt^Yi\i{7^y)dy. (24) 

We note however, that this can also be written as 

1 r°° 
var/~-/ x{CU)^)dx (25) 

TT^ JO 

where C(/)(x) = /o°° /(y) cos{xy)dy denotes the cosine transform of /. This is an exercise 
involving the properties of the Mellin transform, and we leave it to the reader. 

To compute the distribution function, we first turn our attention to the case where 
z/ = —1/2. Our operator B^lcr) has kernel 

— / a {t / a) cos xt cosy tdt 
n JO 

1 r°° 
= — a (t/ a) {cos{{x — y)t) + cos{{x + y)t)) dt 

TX JO 

= '^(C{a){{x~y)a) + C{a){{x + y)a)). 

TT 

This is unitarily equivalent to the operator on ^2(0, a) with kernel 

1 



TT 



.[C{a){x-y) + C{a){x + y)). (26) 



The operator with kernel -(C(cr)(x — y)) is the finite Wiener-Hopf operator, usually 
denoted as Wa{(T), and the operator with kernel -{C{a){x + y)) is the Hankel operator 
Ha{a). (The only difference between this definition of a finite Wiener-Hopf operator and the 
one given earlier for A^ is the difference in the domain. The two are unitarily equivalent.) 
If we consider the operators on ^2(0, 00) in what follows, we will denote them by W{a) and 
H{a) respectively. Also, whenever it is necessary to consider the extension of a to the entire 
real axis, it will always be the even extension. 

7 



Thus the problem of finding the distribution function asymptotically becomes the same 
as computing the Fredholm determinant det(/ + Ba{<j)) = det(/ + Wa{(j) + Ha{(j)) asymp- 
totically. To do this we need some basic facts about Wiener-Hopf operators and we collect 
them in the following theorem. These are well-known and can all be found in [^]. 

Theorem 3 a) Suppose and ip o^^e even bounded functions in Li(R). Then 

W{(I))H{iIj) + H{(I))W{iIj) = H{(f)-^) 

and 

W{(j))W{ip) = W{(f)ip) - H{(f))H{iP). 

b) Suppose (f) and ip are bounded functions in Li(R). // the Fourier transform 0(x) 
vanishes for X negative, then W {iIj)W {(p) = W{(f)ip) and if (j){x) vanishes for x positive, 

thenW{(l))W{'ilj) = W{(l)^). 

We define W{a) and H{a) with a = 1 + / and / in Li by W{a) = I + W{f) and 
H{a) = H{f). Both of these definitions are natural when thought of in a distributional 
setting, and the above theorem holds with these definitions as well. 

The next theorem is of primary importance in the computations that follow. 

Theorem 4 Suppose = 1 + /, (f)~^ = 1 + g where f and g are bounded even functions. 
Then the inverse of W {(j)) + H {(j)) is W{(t)~^) +H{(p-^). 

Proof: Using Theorem 3 parts a) and b) we have, 

= I - H{<P)H{<p-') + i/(00-i) + H{<t>)H{cp-^) 
= I + H{1) = I. 



The same computation holds for {W{(j)^^) + H{(f)~^)){W{(j)) + H{(f))), and so we have shown 
that these operators are inverses of each other. 

It is well known from the theory of Wiener-Hopf operators that under appropriate con- 
ditions det(J + Wa{cr)) has the asymptotic expansion G{(T)°'E{a) where 



G(a) = exp^ riog(l + a(0)rfe 



and E{a) = det {W {(j))W {(p^^)) with = 1 + a. This is simply another version of Theorem 
2. With additional assumptions on 0, it is very easy to adapt this proof to the Bessel case 
z/ = — 1/2 to show that 

det(/ + W^ia) + H^{a)) ~ G{aYE\a) (27) 

and E\a) = det{(W{4>) + H{(j)))W{(f)^^)). Thus to compute the distribution, we need to 
know the form of the above determinant. This is contained in the next theorem. 



Theorem 5 Suppose a = e^^f — 1 where f is even, continuous, piecewise C^ and vanishes 
at infinity. Suppose also that / G Li and the function 

e-(i+a(ir(oi + i/'(op)eL2. 

Then as a -^ oo^ we have 

det{I + W^{a) + H^{a))r^exp{- rkf{x)dx+-f{Q)---J x|C(/)(x)|2rfx}. (28) 

vr JO 4 /vr^ Jo 

Proof: The conditions on a ensure that the above integrals converge, and that the operators 
if (0) and H{f) are trace class. The reader is referred to for details. These assumptions 
also guarantee that (27) holds. It is also easy to see that (?(</)) = exp{^ /Q°°i/c/(a;) dx}. To 
complete the proof we need a concrete representation for (\.ei{{W {(f)) + H {(j)))W {(p^^)) . Define 

h{k) = \ogdet{{W{(f)) + H{(f)))W{(j)-^)) 

where = e*'^-^. Let h{k) = logdet((Vr(0) + if (0))Vr(0_i)). We need to show the second 
derivative of h is constant in k. A standard formula [§ yields 

.'(.) ^ tr((H^(r-))-(H-(0)+^(0))-x -^(^(^) + ^(^))^(^") ) 

x{{w{<p) + H{<p))wir'Hf)) + w{<i)tf)w{r') + H{<ptf)w{r')} 
= tT{{wi<p-')r'w{<p~\-tf)) + iwi<p-')r'w{<p~')wm)wi<p'') 

This uses Theorem 4. Simplifying further and using the fact that H{(j)^^) is trace class we 
have 

h'ik) = tT{iwi(t>-'))-'wi<t>-'Hf)) + wm)wi<t>-') 

Now apply Theorem 3, part a) and the fact that tT{AB) = tT{BA) to find 

-(w{<p-'))-'w{4>-\-zmw{4>-')r'w{r'Hm. 

The conditions on guarantee that the function has a factorization = {g_ + 1 ) {g^ + 1 ) 
such that the Fourier transforms of g^ and g^ vanish for positive and negative real values 
respectively. Then using Theorem 3, part b),it is easy to see that we can write 

W{(P) = W{g. + l)W{g+ + 1), Wir'r' = W{g^ + l)W{g^ + 1). 

A repeated application of these identities allows us to write h"{k) = tiH{if)H{if), and h"{k) 
is independent of k. Thus at this point we have h{k) = ak'^+bk+c where 2a = — tr((ii(/))^. A 
direct computation shows that a = —^ J^ x\C{f){x)\'^ dx. To compute b, notice that h'{0) 
is trii(^/) = 2^ /o°^ C{f{x)) dx. Also h{0) = tr log(i) = 0. Thus the last theorem holds. 



4 The General Case 

In this section we show that under certain conditions, the distribution function for general v 
has the same form as in the case of z/ = — 1/2. The only difference is in the mean which was 
computed in the last section. The attack on the problem is entirely different here. Instead 
of computing determinants asymptotically, we compute the traces of the operators (i?Q,(a))" 
and then piece together the answers to get an answer for the trace of log(/ + Bj^a)) and 
from that to the desired determinant. 

To begin we need to show that ti f{Ba{cr)) makes sense for a class of analytic functions 
/. Just as we can associate the Wiener-Hopf operator with the Fourier transform and a 
multiplication operator, we can also write 

B^{a) = PHM^U 

where H is the Hankel transform and P is the projection on L2(0,l). Since the Hankel 
transform is unitary on ^2(0, 00) (||ll|), the operator norm ||i?cj(o")|| is less than the infinity 
norm ||(t||oo of a. Thus /(i?Q,(cr)) is defined for / analytic on a disk centered at the origin 
with radius ||cr||oo + 6,6 > 0. The operator Bai^cr) is also trace class for a in Li by Mercer's 
Theorem (P Ch.III) as is f{Ba{cr)) for / satisfying the above and /(I) = 0. 

We need some lemmas that will prove to be useful. These may be known already, but 
we include them for completeness. 

Lemma 6 Suppose — l<p<l, 0<A,(5<l,/x<0,p + ;U + (5<0 and < t < 1. Then 

roo 

/ s'P{l + sY\l-s\'^+^\l~ts\-^+^ rfs< Amax(|l-tr^+^,|l-tr^+^)max(r^,rP"^) (29) 
Jo 

where A is some constant independent oft. 
Proof: We have 

roo 

/ sP{l + s)''\l-s\-^+^\l-ts\-^+^ ds 
Jo 

= t-^+^ t s'iX + sY\^ - sl-'+^ll/t - s\-^+^ ds 
Jo 

+t-^^^ f * 5^(1 + sY\i - sr'+^ii/t - s|-i+^ ds 

roo 

+t~^+^ / sP{l + sY\l - s\-'+^\l/t - s\-'+^ ds. 
Jl/t 

We consider each of the above integrals. In each, A is a possibly different constant indepen- 
dent of t but can depend on the other parameters. First, 



s\-'+' ds 




[' sP(l + sY\l-s\-'+^\l/t 
Jo 

Jo 



s\-'+^ds 

< A\t - i\-^+H^-^ . 

10 



Next, 



f * 5^(1 + sY\l - s\-'^+^\l/t - s\-^+^ ds 






Finally, 

/ 5^(1 + SY\1 - s|-^+^|l/t - s\-^+^ ds. 
Jl/t 

< ii-i/tr^+M sp+^^ii/t - s\-^+^ ds 
Ji/t 

__ AU 1 |-l+Ai-p-/^-5-A+l 

Putting this together we have that the original integral is bounded by 

Amax(|l - t|-^+Ml - t\-^+^) max{l, r\rP~^). 

Lemma 7 Suppose —l<p<l, 0<\,6<l,fj,<0,p + fj, + \<0 and t > 1. Then 

roo 

/ sP(l + sVll - s|-^+^|l - ts\-^+^ ds < Amax(\l - t|-^+\ ll - tr^+^) max(l, t'^) (30) 
Jo 

where A is some constant independent oft. 

Proof: The proof of this is almost identical to the previous lemma, and we leave the details 
to the reader. 

Lemma 8 Suppose \x\ < l,Rec > 0,Re{c — b) > 0, and Re{c — a — b) < 0. Then the 
hypergeometric function F{a, b, c, x) satisfies the estimate 

\F{a,b,c,x)\ < A|l-a;|^^(^-"-'') 



with A independent of x. 

Proof: The hypergeometric function satisfies the identity F{a,b,c,x) = (1 — xY~°'~^F{c — 
a,c — b, c, x). Using Euler's integral formula for F, we have 

F(c -a,c-b,c,x) = ^,,.^f^ ,, /' t-''-i(l - tf-\l - tx)-'^+^ dx. (31) 

1 (o)l (c — b) Jo 

11 



The last integral is bounded by J^ t^"" ('-'-'^ {1 - t)R^e{a+6-c-i)^^ ^^ r(Re (c-b))r(Re (a+fe-c)) _ 

J- ( -L VjKj Qi ) 

We next find an integral expression for the trace of (i?Q,(cr))". We proceed informally at 
first and later state things rigorously. Using (|1^) we can write this trace as 



roo rl rl " 



/ • • • / / • • • / n SiXia{xi/a)J„{xiSi)J„{xiSi+i) dsi . . . ds^dxi ...dxn 
JO Jo Jo Jo fj{ 

where si+„ = si. Let a be the Mellin transform of a where c > 0. Then the above becomes 

]^ rc+ioo rc+ioo roo roo rl rl " 

TTT-^ •••/ / •••/ / •••/ \{{six\'^'Ju{xiSi)J^{xiSi+i)a{zi)} 

[2m)"- Jc~ioo Jc-ioo Jo Jo Jo Jo ~~^ 

^^zi+...zn ^^^ ^ ^ ^ dSndXi . . . dXndZi . . . dZn- 

Now use the formula ^ 

x~^Ju{ax)Jy{hx)dx 



-FiuA — ,z/+-:2z/+l: 



2^(a + 6)2--^+ir(l + z/)r(l/2 + |) ' 2 ' 2' '(a + 6)2'' 

where F{a,b]c; z) is the hypergeometric function 2-^1, n times in the integral to get the 



expression 

][ rc+ioo rc+ioo rl rl 



(27rz)" Jc-ioo Jc-ioo JO JO 

X Yl cr(zi 



^z,+...+z„ 



^ _ ^ sf^^V{u + 1 - z,/2)F{u + 1 - 2,/2, u + \-2u + l- j^f^, 



xdsi . . . dsndzi . . . dZn- 
Next we make the change of variables 

51 = s[ 

/ / 

52 — "525]^ 

Sn = ^n • • • ^l 

and the integral becomes 

+ioo rl r— c ^ 

(a/2)"i+-+""2"a(2„) 



1 


J c—ioo 


rc+ioo 
Jc-ioo . 


Jo . 


io 




1 


Jo 


••Sn-1 


{2'n%Y , 





r(z/ + 1 - Zrj2) 

r(i + z/)r(^„/2) 

xs^i+-+^"-i(l + s„ . . . S2)-''^+^"-'F(z/ + 1 - zj2, u+\-2v+l- ^'"■■■'^ ^ 



2' '(l + Sn...S2) 



2' 



'3"(^t)r(z/+ 1 - Zi/2) 2u+l+z,+ i+...z„-i^^ _^ ^-2u+Zi-2 



X IK Vn V ..^w. -;-- .-r^"""^"-^"-Hi + ^.+1 
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xF{u + 1 - Zi/2, z/ + -; 2z/ + 1; -— ^i— -)}(is„ . . . dsidzi . . . dzn- 

2 [l + Si+iY 

Write the inside integral as 

1^ ■ ■ ■ / ^^'^""^ . . . dSn . . .di — / . . . / ... dSn ■ ■ ■ dsi 

Jo Jo Jo Jo Jo 

1 roo roo 



ri poo /"OO 

+ / / . . . / . . .dsn- ■ -dsi. 
Jo Jo Jo 



The last integral in the above sum, inserted in the main integral, is the same as tr (5^(0""))- 
After reversing the order of integration to dsi . . . dSn, the first two terms combine to yield 
limits of integration 



oo roo roo /■! 



"'0 "'0 Jmin(l,^,..., — 1 — ) 

and then the first integration can be done. The result is that 

tr(5,(a))" = tr5,(a") + C(a) 
where C((t) is given by the expression 



(27ri)« }c-ioo • • • Jc-^oo ^ ' ' fj r(l + u)V{Zi/2) 



1 - (min(l, ^, . . . , —^—)Y^+-+^-n 

X 



OO /'OO _ _ 

" S2 ' ' ' ' ' S2...Sn 



Jo Zi + Z2 + . . . + Zn 

X (1 + s„ . . . S2r'''^'--'F{u + 1 - zj2, u+\-2u + Y . ^""•••"^ 

2 (1 + S„...S2 



\2' 



[IT sT'-'^'^^-'-^l + 5.)-^'^+-—^ X F(z. + 1 - ^,„i/2, u+'--2u+l, -4^; 

i=2 2 (1 + SiY 



xdsn ■ ■ ■ ds2dzi . . . dZn- 
We next write this integral as 

]^ rc+ioo rc+ioo foo foo 



rc+ioo rc+ioo roo roo 

/ . . . / G{zi) / . . . / H{zi; Si)dsn . . . ds2dzi . . . dZn- 

Jc—ioo Jc—ioo Jo Jo 



(27ri)" Jc-ioo Jc-ioo Jo Jo 

The idea from here on out is to evaluate this integral asymptotically using complex analysis. 
This will be done in several stages and by breaking the integral into several parts. To begin 
we first consider the interior integration 



roo roo 

/ . . . / H{zi; Si)ds2 ...dsr. 
Jo Jo 



Consider this as an integral over Ri U R2 where i?i is a union of disjoint sets, -Ri = U"^2^j 
such that on Ui, Sj is bounded away from 1 and where R2 is the complement of Ri. 
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Lemma 9 Suppose that —2v — 1 < 0. The integral of H{zi\ Si) over Ui is bounded and the 
Zi variables can be changed in such a way so that the integrated function is analytic in a 
particular z variable to the left of the imaginary axis. 

Proof: For convenience let i = 2 (although the proof is the same for any i) and let 2:1 + ... + 
Zn = z^ with the other variables remaining the same. Suppose that Re 2;, = c for i = 3, . . . , n 
and that c > 0. Suppose also that Rez^ = 6 with |6| < c. We now refer to z'^ as z. Our goal 
is to show that this integral is bounded and that as a function of z is analytic to the left 
of the imaginary axis. By repeated application of Lemma 8, we can say that the integral is 
bounded by a constant times 

/•oo .00 1- (min(l,^, ...,^^))^i+-+^" 

\S2-\\>B Jo Jo Zi + Z2 + . . . + Zn 

n 

i=3 
XS^2'=(1 + S2r^''-'\1 - S2\'-'\l -S2... Snr' 
XdS2 ■ ■ ■ dSn- 

This is valid as long as 2z/ + 1 > and Rezj — 1/2 < 0, which is the case here if we assume 
that c is small enough. Next, we estimate 

1 - (min(l, -,..., ^^))^i+-+^n 

Zi -\- Z2 -\- ■ ■ ■ -\- Zn 

by using the fact that |1 — a;^| < max |z|x | lna;| where the max is taken over the z' values 
on a line connecting and z, x between and 1. Thus, |1 — x^| < K\z\x x~'^ for some 
positive e chosen shortly. Inserting this in the integral we have that the integral is bounded 
by a constant times 

/■ /"OO /"OO ^ 

/ , / .../ E{max(l,(s2...S,r(s2...S,)-^+^)} 

Js2-i>bJo Jo ~L 



i=2 

n 



xY{{i + s,)-'^-\st-'^''\i-sr' 

i=3 
X4-2C(1 + S2)-'''-'\l - S2t'\l -S2... Snr' 

XdS2 ■ . ■ dSn- 

The reason for both terms in the "max" part of the integral is that b could be either positive 
or negative. Now lets begin with the s„ integration. Then the first interior integral has the 
form ^ 

/•oo 
/ S^(l + Sn^-'ll - Snr'\l -82... Snl^^'dSn- 

Jo 
The value for p is either ±a where a = \b — e\ < c. The next step is to apply Lemmas 6 
and 7. We use \ = 5 = c and p as above. The result is that this integral is bounded by a 
constant times 

|1 - S2 • • • Sn-il'^"^ X max(l, (sa . . . S„-l)"^ (S2 • • • Sn-l )"''"'' (S2 • • • Sn-lY^). 

14 



We collect powers and use the lemmas twice with respect to the Sn-i integration and powers 
oi p = ±(2c). At the next integration step the powers oi p = ±3c and so on until we arrive 
at the S2 integration. Here we will have 

Sl\l - S2\'\l + S2\~^''~'dS2 
\S2~1\>B 

where p and q are appropriate powers. These integrals satisfy all the conditions necessary 
for the lemmas as long as c and b are small enough. We will have at most 2" integrals in this 
process. Hence the integral of H over f/j is analytic in the z variable in a strip |Rez| < c by 
the application of Morera's Theorem and Fubini's Theorem. 

We remark here that this proof also is easily modified to show that the interchange of 
integrals done at the beginning of the section are valid and the expression C(cr) is the one 
of interest. 

Lemma 10 Suppose that a has [v\ + 2 derivatives all in Li and that that —2v — 1 < 0. Then 
the integral 

/ . . . / G{zi) ■■■ H{zi; Si)ds2 . . . dsndzi ...dzn 

J c—ioo J c—ioo J J R\ 

is 0{a^^) where 5 > 0. 

Proof: Note that the condition in the hypothesis implies that 

c+oo 

1^(^)112^+^/2^00. (32) 

— ioo 

We first replace the inside integral with a sum of integrals over f/j. For each of these 
we change variables as in the last lemma. We can then perform the integration over the z 
variable by moving it to a line to the left of the imaginary axis. Thus we have that each of 
these integrals is bounded by a constant times 

(7r)"2'' A-ioo Jc-ioo "Jc-ioo il T{l + iy)T{zi/2) ' 

.^{z-j:^^^,Zj)T{u+l-{z-j:z,j^zZj)/2) 

X z— — — — — — — — \dzaz2 . . .azn- 

Til + v)T[{z -Y.z,^z Zj) /2) 

This last integral is bounded by a product of integrals all of the form 

1^(^)11 — T^no\ — 1^^ 

c—ioo i \Z j Lj 



and these in turn are bounded by (|3^) using the basic asymptotics properties of the Gamma 
function. 

We now turn our attention to the region R2. To begin we make another change of vari- 
ables. 

- = 1-4 (33) 

-52 

15 



S2- ■ -Sn 



= l-So-s'o (34) 

S2S3 

: (35) 

1 

-- 1 - 4 - 4 - • • • - 4- (36) 



Under the change of variables, the region R2 is transformed to a region R3 which can be 
assumed to be a symmetric region containing the origin, and where the sum 

\s2 + . . ■ + Sj\ < a < 1 

(we drop the "primes" again) for some a. Notice that the exact form of Ri was unnecessary 
in the previous computation. Thus the integral over R2 is transformed to 

where 



Zl+... + Z„ 



I{Zi\ Si) 



1 - (1 - max(0, g2, • • • , -§2 + ■ ■ ■ + Sn)) 

Zi + ... + Zn 

vie p2-l |„ |2„_i-l| I I \z„-l 

X |S2| ■ ■ ■\Sn\ \S2 + ■ ■ ■ + Sn\ 

y<f{s2,...,Sn,Zi,...,Zn) 

where the function / is smooth in the s variables. 

The following lemmas will help keep track of the contribution of the R^ integral. 

Lemma 11 Suppose Rezi = c,0 < c < 1, for i > 3. Then the integral 

\S2\ I S3 1 ■ ■ ■ \Sn\ \S2 -T ■ ■ ■ Sn\ "S2 • • • USn 

R3 

can be thought of as an analytic function in the z\ variable that can be extended to a strip 
containing the imaginary axis. 

Proof: First note that the following integral with z and w real and between zero and one 
satisfies 

\z-\-w—l 



I Ixl^'^lx + yl'^'^dx < A\y\ 

J a 



where the constant only depends on the z and w variable. A repeated application of this 
estimate in the above integral yields a final integration of 

Thus, once again the analytic continuation argument holds. 
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Lemma 12 Suppose Rezi = c,0 < c < 1, for i >2 and Rezi = d. Then the integral 

\q pile l^2|„ |23-1 \q |2n-l-l|„ _|_ e P"~l/7e rl a 



can he thought of as an analytic function in the zi variable that can be extended to a strip 
containing the imaginary axis. 

Proof: We begin the integration just as in the previous integral After n — 3 integrations we 
arrive at an integral with an estimate of the form 



J a J a 



We can estimate this by looking at three integrals 

fh rl 



a J-l 

b fh/s2 

/ \o l^+("^2)c| ici I 1 |(n-3)c-l 7 7 

/ P2| IS3I |S3 + i| aS3aS2, 

la Jl 

and 

Ja Ja/s2 

We can say, for example, that the last integral is less than a constant times 

and thus is finite for Re zi in a strip about the imaginary axis. The other two integrals are 
handled in the same manner. So by our standard argument the analytic extension is defined. 
Now let us return to our function I{zi; Si). We can write the expression 

1 - (1 - max(0, S2, . . . , S2 + . . . + s„))^i+-^" 



Zi + ...+ Zn 

as 

max(0, S2,...,S2 + ...Sn) + (max(0, S2, • • • , ^2 + • • • s„))^ x g{zi + . . . + z„,S2, ■ ■ ■ Sr 
where the last function is a bounded continuous function in the variables. 
Lemma 13 The contribution of 

]_ rc+ico rc+ioo 



...I G{zi) ... (max(0,S2,...,S2 + •••Sn))^ 

J c—ioo J J Ro 



yZTTlj J c—ioo Jc—ioo J J R2 

X |S2| ■ ■ ■ \Sn\ \S2 -T ■ ■ ■ -T Sn\ 

X/(S2, ■ ■ ■ ,Sn,Zi, . . . , Zn)g{Zi + . . . + Zn, S2, . . . S„) 
X dSn ■ ■ ■ ds2dzi . . . dZn 

IS 0{a~^). 
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Proof: We simply consider the set where say Si + S2 + ■ ■ ■ + Sj is the maximum of the the 
terms. We then expand the square so that we have a term of the form SiSk- We then apply 
the above lemmas after an appropriate re-ordering of the variables and the lemma holds. 

The next step is to replace the function / in the expression for I{zi] Sj) with the first term 
of its Taylor expansion. This expansion gives an "extra" Si (combined with the ones from 
the max(0, S2, ■ ■ ■ ,S2 + ■ ■ ■ + s„)) term in the estimates which, as the above lemmas show, 
is all we need to show that this part of the integral does not contribute in the asymptotic 
expansion. 

So we are finally at the one critical term that gives a contribution in the expansion. This 
term is 

— — — / .../ G{zi) ... max(0,S2,...,S2 + ... + s„) 

X |S2| . . . \Sn\ |S2 -|- • . . -|- S„| 

Xf{0,0,...,0,Zi,...,Zn)dSn. . . dS2dZi . ..dn- 

We can easily compute /(O, 0, . . . , 0, ^i, . . . , Zn) to see that it equals 

" 2-^''-^T{2u + l)T{-Zi/2 + 1/2) 
\^ V{u+l/2)T{u+l-Zi/2) 

We can simplify further using the formula for G{zi) and the duplication formula for the 
Gamma function to arrive at 

^ ' {2-KiY Jc-ioo Jc-ioo ^ ' ' V T{Zi/2) 

X I ... I max(0, S2, . . . , S2 + . . . + Sn) 

X |S2|^^"^ . . . |s„|^"-i"^|s2 + . . . + s„|^"-Ms2 . . . dsndzi . . . dzn + 0(«-^). (37) 

Notice that this expression is now independent of u. Our final steps are to compute the 
contribution from the above integral and we, by the way, finally have an integral which will 
yield a contribution. We begin with a well-known identity due to Mark Kac, which was used 
originally to prove the continuous analogue of the Strong Szego Limit Theorem. It reads 

n 

^max(0, a<^,,a^i +a„^,...,a^^ + ... + a„J = ^^a^^6{a„^ + ... + a^J 

where ^(x) = 1 if a; > and 6{x) = otherwise and the sums are taken over all permutations 
in n variables. 



Because of this identity we can rewrite the integral in (|37D as 



V -±- r^ r^ ial2Y^^-^^-.-'^l^ fr ^(^^)r(-^./2 + 1/2) 
P^ (27ri)" Jc-ioo ■ ■ ■ Jc-ioo ^ ' ' V r(z,/2) 

/ . . . / S2|S2|^^"^ . . . |S„|^"-1"^|S2 + . . . + S„|^"~VS2 . . . dSndZi. . . . dZn- (38) 

J JR3r\{s2+...+Sj>G} 



It is straightforward to see how this identity can be used if the integrand is symmetric in 
the variables. In our case, the integrand is not obviously symmetric in the variables, but can 
always be made so by changing the z variables. Thus we can apply the identity. 

We once again consider the inner integral and call z = Zi + . . . + z„ leaving the other 
variables as is, and show how this inner integral can be thought of as analytic in 2 in a strip 
containing the imaginary axis. The difference is that in this case there will be a pole at 
2 = 0. 

Now we suppose that j > 2. For j = 2 the following computation is almost identical and 
the conclusion is the same. Let us rewrite the inner integral in ( |38D as two integrals 

rb „ 

e Ic \Z — Z2 — Z3 — ...Z„ — 1 I |2n-l — 1|„ _|_ _|_ „ \^n — lj„ J„ 

■52 P2 ■ ■ ■ \Sn\ \S2 ~r ■ ■ ■ ~r Sn\ 0,3^ ■ ■ ■ CIS2 





B 

/e e \Z-Z2-Z3-...Zn-l \ 2„_l-l , , 2n-l^o ^o 

. . . / ■52P2| ■ ■ ■ \Sn\ \S2 -T ■ ■ ■ -T Sn\ CtSn ■ ■ ■ CIS2 

B 

where B is some n — 2 dimensional set. In the first (the computations for the second integral 
being almost identical) of these we make yet another change of variables: 

S3 = S3S2 



^n — ^n^2 

to arrive at 

j ST^ f ... /'|s3|^'"^..|s„|^"^i-^|l + S3 + ... + S„|""-Ms„...rfS3t/S2. 
B/s2 

The original set R^ was chosen to be symmetric and contain the origin. So here we chose 
it to be something convenient, say a cube C with size length /. With this choice we can 
write B / S2 as C / S2 fl {53 + . . . + s„ + 1 > 0}. Next integrate by parts with respect to the S2 
variable. The result is that the above integral becomes: 

slk{s2) - / sldl ds2{kis2))ds2 
Jo 

where 

^(^2) = /... f\s3\'^'K..\Sn\'''-'-^\l + S3 + ... + Sn\'"'^dSn...dS3. 



B/s2 

The function k[s2) has a derivative given by the formula 

k'{s2) = -82^ / /(S3, ...,Sn) {n- S2^{S:i, ..., Sn))dS 

Jd 

where D is the boundary of the set C/s2 which lies in the half-space defined by {53 + . . . + 
s„ + 1 > 0}, the vector n is the outward normal to the surface, the function / is simply the 
one given in the above integral restricted to the surface, and dS is surface measure. We can 
estimate the derivative of k{s2) on any boundary edge to be at most a constant times S2 
for Re Zi = c. Thus we have proved the following: 
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Lemma 14 The function of z defined by 

rb 

■S2P2I '" " . . . |,s„| "'^ P2 + • • ■ + s„| " dSji ■ ■ ■ ds2 

J ' ' ' J 

B 

/•O 

+ / I , I S2\S2\ ■ ■ ■ l-Snl \S2 + ■ ■ ■ + Sn\ CLSn ■ ■ ■ CLS2 

B 

is analytic in a strip containing the imaginary axis except at the point z = 0. Further, the 
contribution of this integral with the z integration moved to a line to the left of the axis is 
given by the residue at z = plus 0{a~^). 

We note here that there are no other poles given our conditions on a , ( ^2]) and the 
formula for G{zi). 

For j > 2, the above computation also shows exactly what the residue is, namely: 

I S3 1 . . . \Sn\ |i + S3 + . . . + s„| asn...ass 

R"-2n{s3+...+Sj>-i} 

I S3 1 . . . \Sn\ I — i -|- S3 -|- . . . -|- s„| aSn . . . as^. 

R"-2n{s3+...+Sj>-i} 

To find an explicit formula for this integral we start with the following formula that can 
be easily proved using formulas for the Beta function. 
For 

< Rep, Reg < 1, Re (p + g) < 1 

\xr'\x + vr^dx = |,|P+.-i2r(p)r(g)cos(vrp/2)cos(vrg/2)_ ^^^^ 

r(p + q) cos((p + g)vr/2) 

Define t{p, q) to be 

2r(p)r(g) cos(7rp/2) cos(7rg/2) 

r(p + q) cos((p + q)Tr/2) 
The residue is then [B is the Beta function) 

n-l j-2 

B{Z2 + . . . + %-l, Zj + ...+ Zn)W t{Zk, Zn + ... + Zk+l) H ^(^fc, Zk+l + • • • + %-l). 

k=j k=2 

We leave this as an exercise to the reader. For j = 2 the residue can also be easily computed 
using the definition of t{p, q) and it is seen to be 

n-l 

11 t{Zki Zn + ■ ■ ■ + ^fe+l)- 
fe=2 

Combining all of the above results we are left with the following theorem. 



oo 
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Theorem 15 Suppose a has [z/] + 2 continuous derivatives in Li. Then 

tr(5„(a)r = tr5„K) + C(a) 
where 

_1 "-1 1 POO 

C{a) = -- y - / xC{a^){x)C{a''-^){x)dx + o(l). 
Proof: Recall we were computing the integral 

V -I^ r^ r^°°ra/2)-+-+-vr-«/2 rr ^(^-)r(V2 + l/2) 

P^ {2'KiY Jc-ioo ■ ■ ■ Jc-ioo ^ ' ' Y T{Zi/2) 

X / . . . / S2|S2|'''"^ • • • |Sn|''""'"^|s2 + • • • + S„|^""^ciS2 • • • dSndZi. . . . dZn- (40) 

J JR3r\{s2+...+Sj>Qi} 

For each j we rename the variables and compute the residue as above. For j > 2 the residue 

is 

-1 f-+i°- f-+i-- :^Jia{z,)T{-Zi/2 + l/2) 



n 



n 



(^27lij Jc-ioo Jc—ioo 2 -'- \^i/^) 

H-Z2 - ■ ■ ■ - Zn)r{{z2 + ... + Zn )/2 + 1/2) 
r((-Z2 - ... - ^n)/2) 



X ;^77 :^-^^ B[Z2 + . . . + Zj-i, Zj + . . . + Zn) 



X n ^(^fc' Zn + ■■■ + Zk+i) n ^(^fc' ^fc+1 + • • • + Zj-i)dz2 . . . (iz„. (41) 

k=j k=2 

Notice that 

t{p,q)t{p + q,r) = 2' 



2r(p)r(g)r(r) cos(p) cos(g) cos(r) 



T{p + q + r) cos((p + q + r)iT/2) 
Using this identity in (§) we have that the above integral is 

„_3 /2 fV a(2;,)r(z,)cos(z,7r/2)r(-2;,/2 + l/2) 



rc+roo rc+ioo "■ 

/ . . . / 2"-37r-"/2 n 

Jc—ioo Jc—ioo ■ n 



(27ri)" ^ Jc-joo Jc-ioo ~^2 r(zj/2) 

-^(-22 - ... - Zn)T{{z2 + . . . + 2;„)/2 + 1/2) 



r((-2;2 - ... - 2;„)/2)r(2;2 + . . . + z^) cos{{z2 + . . . + Zj_i)'k/2) cos{{zj + . . . + 2„)7r/2) 

(42) 



dzo . . . dZm 



From the duplication formula for the Gamma function, this can be simplified to 

1 pc+ioo fC+ioo " 

yZ,nLj J c-ioo Jc— too ■_2 



a{-Z2- ... - Zn){z2 + ... + Zn)sm{{z2 + ... + Zn)n/2) 

X — — — — CiZ2 . . . UiZn. V^^J 

COS[[Z2 + . . . + Zj^i)'k/2) COS[[Zj + . . . + Zn)'K/2) 

Now we change variables with 

Zj-i = Z2 + . . . + Zj-i^ Zn = Zj + . . . + Zn 
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and the above integral becomes 

yZTTtj Jc—ioo Jc—ioo 



]^ rc+ioo rc+ioo J ^ n i 

—-1 ... 2-\-i(n ^{z,))a{z,^, _ . . . _ ^2) n ^(^0 

) Jc—ioo Jc—ioo .•_r, ■_■ 



i=2 i=j 

siniiZj^l + Zn)TT/'2) 
X a{Zn - ... - Zj)a{-Zj_i - Zn){Zj-i + Zn) -. --T -. 77TV"^2 • • • dZn- (44) 

cos(Zj_i7r/2) cos(2;„7r/2) 

The convolution theorem for the Mellin transform shows that this can be reduced to the 
integral 

]_ /-c+ioo rc-\-ioo 



(27ri 



|2 , 

' ' C-lOO J C-lOO 



X (3-(-Zj_i - Zn){Zj_x + ^n) 7 ^"T^T — ^7 j—dZj-idZn. (45) 

cos(^j_i7r/2) cos(z„7r/2) 



Notice this can also be written as 

^c+joo I'C+ioo 



{2m] 



2 , 

' C-lOO J C-lOO 



2-^7r-VJ--2(^^._i)a-^-+i(^„) 



X a{-Zj^i - Zn){zj_i + 2;„) ^ ^ dzj^idzn (46) 

cos(2;j_i7r/2) 



/•C+JOO fC+tOO 

/ / 2-V-Vj-2(^^_i)a"-i+i(2„ 

Jc—ioo Jc—ioo 



X 0-(-2;j_i - Zn>[Zj-i + Znj Z j—dZj^idZn. (47) 



(27ri 

sin(2;„7r/2) 
cos(2;„7r/2) 

Before we proceed further we need three formulas from the theory of Mellin transforms. 
These are _ 

the Mellin transform of / (j){x)dx = z~'^^(z + 1) 

Jx 

where $ is the transform of 0, 

the Mellin transform of X(f)'{x) = —z^{z) 
where $ is the transform of (p, and finally 

2 roo ^ rc+oo 

- xC{(l)){x)C{ilj){x)dx = — :/ ^{z)'^{-z)ztcin{z7T/2)dz. 

Vr JO 2TXi Jc-ioo 

These can be found in any standard table of transforms, although the third requires a 
straightforward computation combined with the convolution theorem. 

So now we apply the second formula along with convolution with respect to the Zn variable 
and we have for each 2 < j < n 

c+ioo —r-, -, ,, , sin(2;,-_i7r/2) , ,, , 

a^~^{zj^i)xa''-^+^a'{-Zj^i) / ' ' dzj^i (48) 



STT^i Jc-ioo cos{zj^in/2) 
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1 /-c+ioo _^ sin(2;,„i7r/2) , , , 

Svr^z Jc-joo cos(2j_i7r/2) 

Next apply the first formula after inserting a factor of Zj^i/zj^i to write the above as 

xC{a^-^){x)C{ a"-^+^a'){x)dx (50) 



27r2 7o 

1 /"CXD /"CXD 

+ 7r^/ xC(a"-^+^)(x)C(/ (T-'-V)(x)t/a; (51) 

27r^ JO A 

or 

+ TTT^ / xC{a'-^){x)C{a^-'+\x) dx. (53) 

ZTT^ J — i JO 

We can do the j = 2,j = n cases separately just as easily (the above formulas are not even 
all required in that case) and putting the two cases together and reindexing when necessary 
we arrive at the conclusion of the theorem. 

Our final step is to extend this to functions other than powers. The standard uniformity 
arguments used in the Wiener-Hopf theory apply here if we can show that 

||tr/(E„(or))-tr5,(/(a))||i = 0(l) 

uniformly for a replaced by 1 — A + Aa and A in some complex neighborhood of [0,1]. 
The details of this are found in [|1^. The norm above is the trace norm. Given sufficient 



analyticity conditions on /, it is only necessary to prove \\Ba{o'i)Ba{cr2)—Ba{ciO'2)\\i = 0(1) 
where the 0(1) here depends on properites of o"j. A trace norm of a product can always be 
estimated by the product of two Hilbert-Schmidt norms and in this case we need to estimate 
the Hilbert Schmidt norm of the operator with kernel 

POO 

-^(ioo)(-2) / ai{t/a)y/xztJ^{xt)Ju{tz)dt. 

JO 

Using integration by parts, and integration formulas for Bessel functions this is easily esti- 



mated to be bounded. For analogous details see [jT4|. Thus for suitably defined / we can 
extend our previous theorem to the more general case. The / of interest is log(l + z). This 
will satisfy the necessary analyticity conditions if we consider small enough k. The necessary 
conditions are collected in the following: 

Theorem 16 Suppose f is a real-valued function with [u] + 2 derivatives all contained in 
Li. Then for sufficiently small k (say k < \\<j\\^) 




ikv ^, , k 
2 27r^ JO 



^{k) ~ exp <! - / ikf{x)dx - -^f{0) -TTY i xC{ff{x)dx 



Proof: The form of the answer follows from the computation of the mean given earlier and 
from the fact that the constant term in the previous theorem is exactly half of the answer 
in Szego's Theorem. Thus the above answer for the log function must be half as well. 

The author would like to thank both Craig Tracy and Harold Widom for many useful 
and helpful conversations. 

23 



References 

[1] E. L. Basor, C. A. Tracy. Variance calculations and the Bessel kernel, J. Statistical 
Physics 73 (1993). 

[2] E. L. Basor, H. Widom. Toeplitz and Wiener-Hopf determinants with piecewise contin- 
uous symbols, J. Functional Analysis 50 (1983) 387-413. 

[3] C. W. J. Beenakker. Universality in the random-matrix theory of quantum transport, 
Phys. Rev. Letts. 70 (1993) 1155-1158. 

[4] A. Bottcher, B. Silbermann. Analysis of Toeplitz Operators, Springer, Berlin, 1990. 

[5] I.e. Gohberg, M.G. Krein. Introduction to the Theory of Linear Nonself adjoint Opera- 
tors, vol 18, Translations of Mathematical Monographs, Amer Math Soc, 1969. 

[6] K. Johannsson. On Fluctuations of Eigenvalues of Random Hermitian Matrices, 
preprint. 

[7] M. Kac,. Toeplitz matrices, translation kernels, and a related problem in probability 
theory. Duke Math.. J. 21 (1954) 501-509. 

[8] M. L. Mehta. Random Matrices, Academic Press, San Diego, 1991. 

[9] P. Sarnak. Arithmetic quantum chaos, preprint. 

[10] A. D. Stone, P. A. Mello, K. A. Muttalib, and J.-L. Pichard, Random theory and max- 
imum entropy models for disordered conductors, in Mesoscopic Phenomena in Solids, 
eds. B. L. Altshuler, P. A. Lee, and R. A. Webb, North-Holland, Amsterdam, 1991. 
(Ch. 9, pp. 369-448.) 

[11] A. Unterberger, J. Unterberger. La Serie discrete de 5*^(2, R) et les operateurs pseudo- 
differentiels sur une demi-droite, Ann. Scient. Ec Norm. Sup. 4 serie, 17, (1984) 83-116. 

[12] C. A. Tracy, H. Widom. Introduction to random matrices, in Proc. 8th Scheveningen 
Conf., Springer Lecture Notes in Physics, 1993. 

[13] C. A. Tracy, H. Widom. Level spacing distributions and the Bessel kernel, Commun. 
Math Phys. 161 (1994) 289-309. 

[14] H. Widom. Szego's limit theorem: the higher-dimensional matrix case, J. of Func. Anal. 
39 (1980) 182-198. 



24 



